System and a method for camera motion analysis and understanding from a video sequence

ABSTRACT

In the event that a moving body (e.g. a person, a car, etc.) is outfitted with a video camera or with a camera-equipped device (e.g. a tablet or a mobile phone), the system described in one aspect is able to understand the motion of the moving by analyzing the video frame sequence captured by the camera. This means that the system can categorize the motion of the body-carrying camera to one of several types (e.g., is this a person walking? is this a person running? etc.), understand the nature of the moving body holding the camera-equipped device (e.g. Is this a car?, Is this a person? etc.) and even to identify the moving body (which car?, which person? etc.).

CROSS REFERENCE TO RELATED APPLICATIONS

This application is a Continuation of U.S. patent application Ser. No. 14/321,424, filed Jul. 1, 2014, now U.S. Pat. No. 9,508,026, entitled “System And A Method For Camera Motion Analysis And Understanding From A Video Sequence,” the entirety of which is incorporated herein by reference.

BACKGROUND

1. Field

An exemplary embodiment of this invention relates to the field of motion analysis and understanding. More specifically one exemplary embodiment relates at least to a method and a system capable of analyzing and understanding the motion that a body, carrying a device equipped with a camera of any kind, undertakes using the video sequence that the camera captures. Another exemplary embodiment of this invention relates to a method and a system capable of distinguishing the environment that the body carrying the camera-equipped device is in and even identifying the person that operates the camera-equipped device.

2. Background

Camera motion analysis and understanding is an important part of understanding video content, and plays a significant role in video browsing, retrieval, editing, printing, etc., in many multimedia systems, including personal computers (PCs), stationary or portable digital entertainment systems, cameras, and mobile devices such as smartphones, tablets, etc.

Existing approaches to motion analysis and content understanding are too slow for common processing systems like PCs and embedded systems like these used in smart cameras and smart mobile appliances (smartphones, tablets, or the like). Existing approaches are typically designed for specific tasks, e.g., tracking the movement of a person (with a known-face model) or a car (with a pre-defined car model), and because of these simplifications have a limited general applicability.

SUMMARY

In the event that a moving body (e.g. a person, a car, etc.) is outfitted with a video camera or with a camera-equipped device (e.g. a tablet or a mobile phone), the system described in one exemplary embodiment of the current invention is able to understand the motion of the moving by analyzing the video frame sequence captured by the camera.

This means that the system can categorize the motion of the body-carrying camera to one of several types (e.g., is this a person walking? is this a person running? etc.), understand the nature of the moving body holding the camera-equipped device (e.g. Is this a car?, Is this a person? etc.) and even to identify the moving body (which car?, which person? etc.).

Furthermore, in the event that a person operates the lens of a camera in some way (e.g. the camera operator zooms-in or zooms-out), the system described in one embodiment of the current invention is able to understand the nature of the operator's control command by analyzing the video frame sequence captured by the camera.

In one aspect of the current invention, the camera motion is analyzed through the calculation of camera motion parameters. These parameters are able to describe the motion that an image produced by a specific moving camera undertakes, due to this motion.

The method of transforming camera motion parameters into camera motion information goes through the following exemplary stages:

-   -   Calculation of these motion parameters by: determining local         motion vectors from pairs of successive video frames captured by         the camera.     -   Removing local motion vector outliers.     -   Calculation of the global motion of the camera using these local         motion vectors.     -   Modeling the time-behavior of each global motion parameter using         a set of polynomial time-series models. These models are trained         to model the time behavior of each parameter over predetermined         time segments or continuously in a recursive way, updating         continuously a state to model the current motion.     -   Determination of the motion status from the camera motion         parameters behavior over time using the coefficients of the         polynomial time-series models. Motion status is then obtained by         performing pattern classification and recognition using         supervised learning techniques. Due to the potentially large         dimensionality of cluster samples, dimensionality reduction         techniques (using Multiclass Linear Discriminant Analysis as we         describe herein) are applied in this stage, to dramatically         reduce the computational complexity and make the system         appropriate for embedded system implementation and other         applications where lower computational complexity is beneficial.     -   Categorizing the motion to one or more of various discrete         motion categories and outputting the motion state in the form of         a membership percentage to the various motion categories.

BRIEF DESCRIPTION OF THE DRAWINGS

The exemplary embodiments of the invention will be described in detail, with reference to the following figures, wherein:

FIG. 1 illustrates an exemplary motion analysis system;

FIG. 2 illustrates an exemplary camera motion understanding device;

FIG. 3 illustrates a 3-D space showing the full motion of a camera in the 3-D space characterized by a total of eight parameters;

FIG. 4 illustrates a specific time dependence of the value of the parameter;

FIG. 5 illustrates model parameter vector estimation for each parameter;

FIG. 6 illustrates an exemplary 3-dimensional feature space; and

FIG. 7 illustrates a visualization showing the membership percentages of each motion vector to the various classes.

DETAILED DESCRIPTION

In accordance with one exemplary embodiment, a system is disclosed that is able to understand the motion of a device equipped with a digital camera by analyzing the video sequence captured by this camera, using a system like the one shown in FIG. 1. This system uses a camera sensor (11 in FIG. 1) which captures the video frames (12 in FIG. 1), stores the most recent frame in a memory (13 in FIG. 1) and then processes it with a Motion Understanding device (14 in FIG. 1), comprised of a storage section (15 in FIG. 1) and a processing section (which includes one or more processors) (16 in FIG. 1) and in which the motion state is determined.

An exemplary Motion Detection Device through which the system extracts and manages the optical flow in order to calculate the motion information is shown in FIG. 2.

This exemplary system functions as follows: First two consecutive frames I₁ and I_(i+1), (12 in FIG. 1) are input into the Image Data Input Unit (221 in FIG. 2) from the Storage Memory (13 in FIG. 1) and are temporarily stored into the Input Image Data memory (21 in FIG. 2). The data are then fed into the Global Motion Estimation Unit (222 in FIG. 2). The global motion Estimation unit estimates the global motion of frame I_(i) in relation to the frame T_(i+1), in the video sequence, produces a global motion parameter set corresponding to the physical motion of the frame and stores the data to the global motion data memory (22 in FIG. 2). These data are then fed into the motion-model estimation unit (223 in FIG. 2). This unit estimates a motion model from the time-series corresponding to the development of each parameter over time and stores the models to the Motion-Model Data memory (23 in FIG. 2). These model data are then fed into the Classification unit (224 in FIG. 2).

The Classification unit (224 in FIG. 2) uses classification information stored into classification data memory (24 in FIG. 2) at an off-line training phase and classifies the model data into one of the pre-learned states. Then the data are fed into the Motion-State Output unit (225 in FIG. 2), which in turn outputs the recognized motion state.

In the following sections, the above-referenced units are explained in detail.

Global Motion Estimation Unit (222 in FIG. 2)

The Global Estimation unit is responsible for estimating the motion parameters of a camera from the analysis of two subsequent frames.

The aim of the Global Motion Estimation unit is to analyze optical flow as represented by the “cleared” (i.e. outlier-free) local motion vectors, in order to obtain the motion of the camera.

The full motion of a camera in the 3-D space can be characterized by a total of eight parameters, known as the Degrees of Freedom (DoF). These are two translational components, a rotational component, two scale components, two shearing and a non-linearity component of the shearing. However, in one approach the motion of the camera is estimated using the six most dominant parameters. That is, two Translational components (T_(x), T_(y), 71 in FIG. 3), a Rotational component (θ, 72 in FIG. 3), two scale components (sα in x-dimension and sb in y-dimension, 73 in FIG. 3) and a shearing component (h, 74 in FIG. 3).

A system appropriate to fulfill this task has been described in U.S. patent application Ser. No. 13/952,894 entitled “A SYSTEM AND A METHOD FOR MOTION ESTIMATION BASED ON A SERIES OF 2D IMAGES,” filed Jul. 29, 2013, the contents of which are incorporated herein by reference in their entirety. The output of this unit is a global motion parameter vector PAR={T_(x), T_(y), θ, sα, sb, h} which is stored in Global Motion Data memory (22 in FIG. 2).

Motion-Model Estimation Unit (223 in FIG. 2)

One function of this unit is to estimate a model of the motion of the camera over time and works as follows:

The motion of the camera produces a different set of parameters PAR_(i)={T_(xi), T_(yi), θ_(i), sα_(i), sb_(i), h_(i)} at every time t for each pair of frames I_(i) and I_(i+1). This causes a specific time dependence of the value of the parameter as shown in FIG. 4 as an example. For a certain time period, the motion of the camera-equipped body, corresponding to a specific action (shaking, zooming, walking, driving etc.) or a specific environment (car, train etc.), results in a specific time behavior, which can be modeled by using polynomial modeling.

The models are learned, potentially off-line, as follows:

Initially, a video is recorded, when the camera-equipped device undertakes motion of a specific type (e.g. walking, running, etc.) or the device operates in a specific environment (e.g. home, car, etc.) and/or is operated by a specific person. Then, the global motion estimation block (222 in FIG. 2) is employed to extract the time-series behavior of the camera motion parameters, corresponding to various device motions as described above.

Subsequently each time-series sequence is partitioned to a specific number of time-segments N_(S) corresponding to time intervals of a specific duration. For each time interval, a polynomial time-series model is trained, each model corresponding to a vector. Given the fact that in each of the segments the motion of the camera-equipped device is of the same type, then these N_(S) vectors of model parameters should form a cluster in the model feature space, that is they should form a multi-dimensional Euclidean space for which the model vector parameters are the coordinates.

Repeating this procedure for various device motions, operating environments and users, a separate cluster in the feature space can be created for each. For a number of NM different motions, the result will be the formation of N_(M) clusters, each one corresponding to a different kind of motion.

In the recognition phase, the following two exemplary schemes can be followed:

According to the first scheme the time-series sequence of a parameter can be partitioned to a specific number of segments N_(S) corresponding to time intervals of an appropriate duration. Then, for each time interval, a model is trained and then classified by using a classification scheme described elsewhere, e.g. as in 224 in FIG. 2 as described below.

In the second scheme, which is followed in the current embodiment, a model is trained continuously and recursively: After a first training period, necessary for the system to stabilize its functional state, a polynomial time-series model is produced and updated for every sample of the parameter under consideration which is then classified. This produces a motion state characterization for every sample (parameter value), resulting in a continuous and real-time motion-state estimation system.

Polynomial Modeling

There are a number of methods for modeling the time behavior of a parameter. One of the most common approaches is the polynomial time-series modeling. A polynomial model uses a generalized notion of transfer functions to express the relationship between the input, u(t), the outputy(t), and a white noise source e(t) using the equation [1]:

$\begin{matrix} {{{A\left( {q,\theta} \right)}{\hat{y}\left( {t❘\omega} \right)}} = {{\frac{B\left( {q,\omega} \right)}{F\left( {q,\omega} \right)}{u\left( {t - d} \right)}} + {\frac{C\left( {q,\omega} \right)}{D\left( {q,\omega} \right)}{e(t)}}}} & (1) \end{matrix}$ where, A(q, ω)=1+a ₁ q ⁻¹ + . . . +a _(NA) q ^(−NA) B(q, ω)=1+b ₁ q ⁻¹+ . . . +b _(NB) q ^(−NB) C(q,ω)=1+c ₁ q ⁻¹ + . . . +c _(NC) q ^(−NC) D(q,ω)=1+d ₁ q ⁻¹ + . . . +d _(ND) q ^(−ND) F(q,ω)=1+f ₁ q ⁻¹ + . . . +f _(NF) q ^(−NF)

The functions A, B, C, D, and F are polynomials of various orders expressed using the time-shift operator q, and ω is a vector of the coefficients of the polynomial. In practice, not all the polynomials are simultaneously active; by selecting proper values for the polynomials A, B, C, D, F, simpler forms are employed, such as ARX, ARMAX, Output-Error, and Box-Jenkins models[1].

In various embodiments, any of the aforementioned methods can be used.

The general polynomial equation (1) is written in terms of the time-shift operator q⁻¹. To understand this time-shift operator, consider the following discrete-time difference equation: y(t)+a ₁ y(t−T)+a ₂ y(t−2T)=b ₁ u(t−T)+b ₂ u(t−2T)  (2a) where y(t) is the output, u(t) is the input, and T is the sampling interval. q⁻¹ is a time-shift operator that compactly represents such difference equations using q⁻¹u(t)=u(t−T).

Then (2a) can be written as y(t)+a ₁ q ⁻¹ y(t)+a ₂ q ⁻² y(t)=b ₁ q ⁻¹ u(t)+b ₂ q ⁻² u(t)  (2b)

The model (1) is uniquely defined by the parameter vector ω, which is defined as: ω=[a ₁ ,a ₂ , . . . ,a _(NA) ,b ₁ ,b ₂ , . . . , b _(NB) ,c ₁ ,c ₂ , . . . , c _(NC) ,d ₁ ,d ₂ , . . . d _(ND)]  (3) where NA, NB, NC, ND are the order of the polynomials A, B, C, D.

In a preferred embodiment, the ARMA model is used. An ARMA model can be described by the following equation:

$\begin{matrix} {{y\left( {t❘\omega} \right)} = {{\frac{C\left( {q,\omega} \right)}{A\left( {q,\omega} \right)}{e(t)}} = {w\left( {t,\omega} \right)}}} & (4) \end{matrix}$ if we define φ(t, ω) =[e(t−1) . . . , e(t−NC)−w(t−1,ω) . . . w(t−NA, ω)] then the model (4) can be written as: y(t|ω)=φ^(T)(t, ω)ω  (5)

The goal of using such a model is to have a system that can learn to predict values of the variable y at time t by using a number N of past values of y e.g. y(t−1), y(t−2) , . . . y(t−N). In this invention the variable y is made equal to one of the parameters of the global motion vector PAR which is the output of the Global Motion Estimation unit (222 in FIG. 2) and which is stored in Global Motion Data memory (22 in FIG. 2).

Model Estimation is equivalent to the calculation of the model parameter vector ω. This model parameter vector is computed by requesting the minimization of the estimation error, which is quantified by using a specific error function E(ω)) between the training data sequence and the estimated data. To this end, minimization of the error corresponds to the minimization of the error function: E(ω,N)=Σ_(i=1) ^(N)λ(ε(t,ω)  (6) where ε(t,ω)=y(t)−{circumflex over (y)}(t|ω)  (7) and λ(.) is an error function.

In one implementation the Least Mean Square estimator, λ(ε)=(½)ε² can be used. The error function can also be time depended of the form λ(ε,t). This is useful in cases where various measurements have different reliability, and therefore they can be assigned different weights, or when the model is learned in the progress of time. The time variance can also be incorporated by using a multiplication time-varying function β(N,t). Incorporating this function in (6) and using (7) we get: E(ω,N)=(½)Σ_(k=1) ^(N)β(N,t)[y(k)−φ^(T)(k,ω)ω]²  (8)

By using this form of the error function, the requested parameter set ω can be obtained as {circumflex over (ω)}_(N) ^(LMS) by the following relation. {circumflex over (ω)} _(N) ^(LMS)=argmin_(ω)(E(ω,N))  (9)

Equation (9) means that the requested solution will be calculated as the minimum solution.

Equation (9) can be solved in closed form as follows: {circumflex over (ω)} _(N) ^(LMS) =R ⁻¹(t)fn(t)  (10a) where R ⁻¹(t)=Σ_(k=1) ^(N)β(t,k)φ(k,ω)φ^(T)(k,ω)  (10b) fn(t)=Σ_(k=1) ^(N)β(t,k)φ(k)y(k)  (10c)

In many cases it is useful to have a model of the system available on-line, while the system is in operation. The model should then be based on observations up to the current time and be gradually built in progress of time. The methods for computing online models are called recursive system identification methods and employ an adaptation scheme based upon the on-line data.

In an exemplary embodiment and in a recursive formulation the parameter set at time t can be calculated using the following formulas: {circumflex over (ω)} _(N) ^(LMS)(t)={circumflex over (ω)} _(N) ^(LMS)(t−1)+ R ⁻¹(t)φ(t)[y(k)−φ^(T)(k,ω){circumflex over (ω)}_(N) ^(LMS)(t−1)]  (11a) R (t)=λ(t) R (t−1)+ϕ(t)ϕ^(T)(t)  (11b) where β(t,k)=λ(t)β(t−1,k)  (11C)

The factor λ(t) is in this case an adaptation gain, regulating the rate at which the parameter set {circumflex over (ω)} is adapted in progress of time: Small λ(t) corresponds to slow adaptation and large λ(t) corresponds to a fast adaptation.

The model parameter vector

(43 in FIG. 5) is estimated for each parameter P_(i) (41 in FIG. 5) of the motion vector (41 in FIG. 5) extracted by the motion estimation unit This procedure (45 in FIG. 5) is executed in a recursive way for a specific amount of time, until a steady solution is reached, and after that a set of N_(P) vectors is formed of size MO each, MO being the order of the polynomial model. Finally, a feature vector

={ω₁₁, ω₁₂, . . . , ω_(1M0), . . . , ω_(NP1), ω_(NP2), . . . , ω_(2NPMO)} is formed (44 in FIG. 5), which is stored in the Motion-Model Data Memory (23 in FIG. 2) to be later used by the classification unit (224 in FIG. 2).

Classification Unit (224 in FIG. 2)

This unit is focused on classification of a segment of the time series signal corresponding to the each motion parameter.

In one exemplary embodiment of the current implementation, a pattern classification scheme is used for classification. To this end, the system has been previously trained offline, using a database with models corresponding to specific camera or lens motions (e.g. walking, zooming etc), moving bodies (e.g. car, person, a train etc), or users. For each one, the models described in the previous section are evaluated. The various models (serving the role of “features”) are then combined into a total feature vector. This feature vector is then projected in an Euclidean space (referred as a the “feature space”). This Euclidean space is defined as a multi-dimensional space with as many dimensions as the feature vector. In such a projection, the feature vectors corresponding to specific camera motions are concentrated (clustered) in separate areas of the multi-dimensional feature space. Consider the example shown in FIG. 6 incorporating a 3-dimensional feature space (51 in FIG. 6). Each point in this space is defined by a vector

₁={f₁, f₂, f₃} where f₁, f₂, f₃ are its three coordinates. The projection of the several feature vectors on this axis-system creates two clusters, one corresponding to a motion of e.g type-1 (53 in FIG. 6) and one corresponding to motion of e.g type-2 (52 in FIG. 6).

The next step is to define the centers of the individual clusters. In one implementation this is achieved via the calculation of the center of mass of each cluster. The center of mass has coordinates Ĉ={f₁, f₂, . . . , f_(D)} where D is the dimensionality of the feature space, and each coordinate f _(k) is defined as:

$\begin{matrix} {{\hat{f}}_{k} = {\frac{1}{N_{s}}{\sum_{i}f_{ki}}}} & (12) \end{matrix}$ where N_(S) is the number of samples (regions) participating in each cluster. In the 3-dimensional example referred before, the centers of the clusters are indicated as C1 (56 in FIGS. 6) and C2 (57 in FIG. 6).

When a new sample is tested, its feature vector FV _(T) is obtained from the Motion-Model Estimation Unit (223 in FIG. 2). This corresponds to a point in the feature space. In order to test into which cluster this test point belongs, the distance of this point from the centers of the clusters is computed using some distance measure such as the L1 (or absolute) distance, or L2 (or root mean squared) distance.

In one implementation the L2 distance is used which is defined as follows: in Cartesian coordinates, if P=(p₁, p₂, . . . , p_(n)) and Q=(q₁, q₂, . . . , q_(n)) are two points in Euclidean n-space, then the L2 or Euclidean distance from P to Q, or from Q, to P is given by the following expression: d({circumflex over (P)}, {circumflex over (Q)})=d({circumflex over (Q)}, {circumflex over (P)})=√{square root over (Σ_(i−1) ^(n)(q _(i) =p _(i))²)}  (13)

In the 3-dimensional example of FIG. 5, the distance of the test point T (55 in FIG. 5) from the cluster-center C1 52 in FIG. 6) is d1 (58 in FIG. 6) and from the cluster-center C2 (57 in FIG. 6) is d2 (54 in FIG. 6).

In a different implementation the samples of each cluster can be modeled as multi-dimensional normal distributions N(μ, Σ) having a mean vector μ and a covariance matrix Σ. After doing this, distribution distance measures can be used such as the Mahalanobis distance, Kolmogorov-Smirnov distance the Kullback-Leibler divergence, χ² statistics distance etc[2], in order to calculate the distance of a sample (or a cluster of samples forming a distribution) from a specific cluster.

Once the distances of the test point from the centers of the clusters (or the clusters considered as distributions) are computed, the decision about into which cluster this sample belongs to, is taken according a proximity criterion. That is, the point belongs to the nearest cluster according to the distance measure used. Once this decision has been made, the segment under test has been classified.

In a different implementation and if the dimensionality of the feature space (corresponding to the size of the feature vector) is large, dimensionality reduction techniques like PCA (Principal Component Analysis) or LDA (Linear Discriminant Analysis)[3] or a combination of these two can be used. In this way, the dimensionality of the feature space can be reduced dramatically to a number of N_(M)−1, where N_(M) is the number of different motions that the system will be able to recognize. This fact is extremely important in embedded system implementations where the resources, in terms of memory and processing power, are often limited. In such platforms the linear algebra operations are computationally intensive when engage large matrices. Therefore, reducing the size of the related matrices is a very critical step towards relaxing the computational needs and achieve real-time performance in embedded systems.

In the current embodiment, a Multiclass Linear Discriminant Analysis (MLDA) is used as dimensionality reduction scheme. Multiclass Linear Discriminant Analysis aims to map a set of samples from N_(M) classes into the linear subspace that best separates the samples of the different classes, while keeping each class as compact as possible. In other words, LDA seeks for the linear projection that maximizes the ratio of between-class variance to the within-class variance. The within-class covariance is the covariance of the samples participating to each class and indicates the spread of the class in the feature space having the meaning of the size or the volume of the class. The between class covariance is computed for the total population and indicates the spread of the total population in the feature space.

In the multivariate case, the separation of the classes along a direction w of the feature space is defined as:

$\begin{matrix} {S = \frac{{\overset{\_}{w}}^{T}{\sum_{btw}\overset{\_}{w}}}{{\overset{\_}{w}}^{T}{\sum\overset{\_}{w}}}} & (14) \end{matrix}$

where Σ is the covariance matrix of the dataset, Σ_(btw) is the between-class covariance matrix:

$\begin{matrix} {\sum_{btw}{= {\frac{1}{N_{M}}{\sum_{i = 1}^{N_{M}}{\left( {{\hat{\mu}}_{i} - \hat{\mu}} \right)\left( {{\hat{\mu}}_{i} - \hat{\mu}} \right)^{T}}}}}} & (15) \end{matrix}$ {circumflex over (μ)} is the mean of the total population and {circumflex over (μ)}_(i) is the mean of the i-th class. Seeking of the k∈{1, . . . , N_(M)−1} orthogonal directions which offer the highest class separation is equivalent of solving the following generalized eigenvalues problem: Σ_(btw) Ŵ=Σ _(wn) Ŵλ,   (16) where Σ_(wn) is the average within-class covariance matrix. In order to simplify computations, instead of using (15) the between-class covariance Σ_(btw) is computed by subtracting the within-class covariance from the covariance matrix of the data i.e. Σ_(btw)=Σ−Σ_(wn)  (17)

The k requested orthogonal directions are calculated by selecting the k column vectors of W (eigenvectors), which corresponds to the largest k values from the diagonal of the eigenvalues matrix λ. These form a projection matrix M_(P) of dimension nx(N_(M)−1). To this end, projecting each feature vector {circumflex over (P)}=(p₁, p₂, . . . , p_(n)) to the new subspace is equivalent to a vector multiplication of the vector P with the projection matrix M_(P) resulting in a new vector {circumflex over (P)}′ with reduced dimension: {circumflex over (P)}′={circumflex over (P)}

  (18)

The classification process as described above can then be applied to the reduced-dimensionality feature space.

In a further different implementation, an Artificial Neural Network classifier or any other type of classifier can be used, either on the original or at the reduced feature space.

In an even further different implementation, apart for being able to discretely categorize the motion to various discrete motion categories, the system is able to output the motion state in the form of a membership percentage to the various motion categories. This is achieved by assigning to each sample a score vector

={sc₁, sc₂, . . . , sc_(D)} formed as follows:

First a distance d_(i) of a sample FV_(T) from each class-center Ĉ_(i) is computed, using a multivariate Gaussian function:

d i = 1 ( 2 ⁢ π ) D ⁢  ∑ wn i  ⁢ exp ⁡ ( - 1 2 ⁢ ( T - μ ^ i ) T ⁢ ∑ wn i - 1 ⁢ ( T - μ ^ i ) ) ( 19 ) where λ_(i) is the class mean vector, Σ_(wni) is the within class covariance matrix of each class and |Σ_(wni)| its determinant. Then the elements of the score vector SC are computed as follows: sc _(i) =d _(i)/Σ_(k) ^(D) d _(k)  (20)

The vector elements of the vector

corresponding to the membership percentages of each motion vector to the various classes. A possible visualization exploiting the system output in this case, is shown in FIG. 7. In this, each element of the vector SC, corresponds to a specific bar.

Motion State Output unit (226 in FIG. 2)

One aim of this unit is to output the motion state or the motion state vector to the system output.

The exemplary systems and methods of this disclosure have been described in relation to camera motion analysis. However, to avoid unnecessarily obscuring the present disclosure, the preceding description omits a number of known structures and devices. This omission is not to be construed as a limitation of the scopes of the claims. Specific details are set forth to provide an understanding of the present disclosure. It should however be appreciated that the present disclosure may be practiced in a variety of ways beyond the specific detail set forth herein.

Furthermore, while the exemplary aspects, embodiments, options, and/or configurations illustrated herein show the various components of the system collocated, certain components of the system can be located remotely, at distant portions of a distributed network, such as a LAN and/or the Internet, or within a dedicated system. Thus, it should be appreciated, that the components of the system can be combined in to one or more devices, such as a Personal Computer (PC), laptop, netbook, smart phone, Personal Digital Assistant (PDA), tablet, etc., or collocated on a particular node of a distributed network, such as an analog and/or digital telecommunications network, a packet-switch network, or a circuit-switched network. It will be appreciated from the preceding description, and for reasons of computational efficiency, that the components of the system can be arranged at any location within a distributed network of components without affecting the operation of the system. Similarly, one or more functional portions of the system could be distributed between a camera device(s) and an associated computing device(s).

Furthermore, it should be appreciated that the various links connecting the elements can be wired or wireless links, or any combination thereof, or any other known or later developed element(s) that is capable of supplying and/or communicating data to and from the connected elements. These wired or wireless links can also be secure links and may be capable of communicating encrypted information. Transmission media used as links, for example, can be any suitable carrier for electrical signals, including coaxial cables, copper wire and/or fiber optics, and may take the form of acoustic or light waves, such as those generated during radio-wave and infra-red data communications.

Also, while the flowcharts and methodology have been discussed and illustrated in relation to a particular sequence of events, it should be appreciated that changes, additions, and omissions to this sequence can occur without materially affecting the operation of the disclosed embodiments, configuration, and aspects.

A number of variations and modifications of the disclosure can be used. It would be possible to provide and/or claim some features of the disclosure without providing others.

Optionally, the systems and methods of this disclosure can be implemented in conjunction with a special purpose computer, a programmed microprocessor or microcontroller and peripheral integrated circuit element(s), an ASIC or other integrated circuit, a digital signal processor, a hard-wired electronic or logic circuit such as discrete element circuit, a programmable logic device or gate array such as PLD, PLA, FPGA, PAL, special purpose computer, any comparable means, or the like. In general, any device(s) or means capable of implementing the methodology illustrated herein can be used to implement the various aspects of this disclosure. Exemplary hardware that can be used for the disclosed embodiments, configurations and aspects includes computers, handheld devices, telephones (e.g., cellular, Internet enabled, digital, analog, hybrids, and others), and other hardware known in the art Some of these devices include processors (e.g., a single or multiple microprocessors), memory, nonvolatile storage, input devices, and output devices. Furthermore, alternative software implementations including, but not limited to, distributed processing or component/object distributed processing, parallel processing, or virtual machine processing can also be constructed to implement the methods described herein.

In yet another embodiment, the disclosed methods may be readily implemented in conjunction with software using object or object-oriented software development environments that provide portable source code that can be used on a variety of computer or workstation platforms. Alternatively, the disclosed system may be implemented partially or fully in hardware using standard logic circuits or VLSI design. Whether software or hardware is used to implement the systems in accordance with this disclosure is dependent on the speed and/or efficiency requirements of the system, the particular function, and the particular software or hardware systems or microprocessor or microcomputer systems being utilized.

In yet another embodiment, the disclosed methods may be partially implemented in software that can be stored on a storage medium, executed on programmed general-purpose computer with the cooperation of a controller and memory, a special purpose computer, a microprocessor, or the like. In these instances, the systems and methods of this disclosure can be implemented as program embedded on personal computer such as an applet, JAVA® or CGI script, as a resource residing on a server or computer workstation, as a routine embedded in a dedicated measurement system, system component, or the like. The system can also be implemented by physically incorporating the system and/or method into a software and/or hardware system.

Examples of the processors as described herein may include, but are not limited to, at least one of Qualcomm® Snapdragon® 800 and 801, Qualcomm® Snapdragon® 610 and 615 with 4G LTE Integration and 64-bit computing, Apple® A7 processor with 64-bit architecture, Apple® M7 motion coprocessors, Samsung® Exynos® series, the Intel® Core™ family of processors, the Intel® Xeon® family of processors, the Intel® Atom™ family of processors, the Intel Itanium® family of processors, Intel® Core® i5-4670K and i7-4770K 22nm Haswell, Intel® Core® i5-3570K 22nm Ivy Bridge, the AMD® FX™ family of processors, AMD® FX-4300, FX-6300, and FX-8350 32nm Vishera, AMD® Kaveri processors, ARM® Cortex™-M processors, ARM® Cortex-A and ARM926EJ-S™ processors, other industry-equivalent processors, and may perform computational functions using any known or future-developed standard, instruction set, libraries, and/or architecture.

Although the present disclosure describes components and functions implemented in the aspects, embodiments, and/or configurations with reference to particular standards and protocols, the aspects, embodiments, and/or configurations are not limited to such standards and protocols. Other similar standards and protocols not mentioned herein are in existence and are considered to be included in the present disclosure. Moreover, the standards and protocols mentioned herein and other similar standards and protocols not mentioned herein are periodically superseded by faster or more effective equivalents having essentially the same functions. Such replacement standards and protocols having the same functions are considered equivalents included in the present disclosure.

The present disclosure, in various aspects, embodiments, and/or configurations, includes components, methods, processes, systems and/or apparatus substantially as depicted and described herein, including various aspects, embodiments, configurations embodiments, subcombinations, and/or subsets thereof. Those of skill in the art will understand how to make and use the disclosed aspects, embodiments, and/or configurations after understanding the present disclosure. The present disclosure, in various aspects, embodiments, and/or configurations, includes providing devices and processes in the absence of items not depicted and/or described herein or in various aspects, embodiments, and/or configurations hereof, including in the absence of such items as may have been used in previous devices or processes, e.g., for improving performance, achieving ease and\ or reducing cost of implementation.

The foregoing discussion has been presented for purposes of illustration and description. The foregoing is not intended to limit the disclosure to the form or forms disclosed herein. In the foregoing Detailed Description for example, various features of the disclosure are grouped together in one or more aspects, embodiments, and/or configurations for the purpose of streamlining the disclosure. The features of the aspects, embodiments, and/or configurations of the disclosure may be combined in alternate aspects, embodiments, and/or configurations other than those discussed above. This method of disclosure is not to be interpreted as reflecting an intention that the claims require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive aspects lie in less than all features of a single foregoing disclosed aspect, embodiment, and/or configuration. Thus, the following claims are hereby incorporated into this Detailed Description, with each claim standing on its own as a separate preferred embodiment of the disclosure.

Moreover, though the description has included description of one or more aspects, embodiments, and/or configurations and certain variations and modifications, other variations, combinations, and modifications are within the scope of the disclosure, e.g., as may be within the skill and knowledge of those in the art, after understanding the present disclosure. It is intended to obtain rights which include alternative aspects, embodiments, and/or configurations to the extent permitted, including alternate, interchangeable and/or equivalent structures, functions, ranges or steps to those claimed, whether or not such alternate, interchangeable and/or equivalent structures, functions, ranges or steps are disclosed herein, and without intending to publicly dedicate any patentable subject matter.

References: (All of Which are Incorporated Herein by Reference in Their Entirety)

-   -   1. System Identification: Theory far the User, Second Edition,         by Lennart Ljung, Chapter 4, pages 87-88, Prentice Hall PTR,         1999.     -   2. Empirical Evaluation of Dissimilarity Measures for Color and         Texture, J. Puzicha, et all. The Proceedings of the Seventh IEEE         International Conference on Computer Vision, 1999. (Volume 2),         pp. 1165-1172.     -   3. Pattern Classification, 2nd Edition, Richard O. Duda,         Peter E. Hart, David G. Stork Chapter 4, pp. 44-51 and Chapter         10, pp. 53, Willey 2000. 

The invention claimed is:
 1. A method of classifying a motion status of a moving body comprising: capturing video frames using a camera sensor that is attached to the moving body; determining the motion status of the moving body by determining, by a camera motion understanding device, local motion vectors from pairs of successive video frames classifying, by a classification unit, the motion status of the moving body as being one of a discrete number of motion status options by using mathematical autoregressive modeling of global motion parameters from the pairs of successive video frames, wherein the mathematical autoregressive models utilize a least a mean square estimator (LMS) to recursively determine the mathematical autoregressive model parameters, and wherein the mathematical autoregressive modeling utilizes inputs, outputs and a white noise source, and wherein the mathematical autoregressive model is updated for every sample of the model parameter under consideration and produces a motion state classification for each model parameter by minimizing an error function.
 2. The method of claim 1, wherein the model parameters can be continuously updated over time to reflect the current status of motion of the camera.
 3. The method of claim 2, wherein the model parameters are each of: two translational components, one rotational component, two scale components, two shearing and one non-linearity component of shearing.
 4. The method of claim 2, wherein the model parameters are each of: two translational components, one rotational component, two scale components, and one shearing component.
 5. The method of claim 4, wherein the error function is time-dependent.
 6. The method of claim 5, wherein the error function associated with the LMS estimator is further adjusted for being responsive to various time-dependent conditions.
 7. A system for classifying a motion status of a moving body, the system comprising: a camera sensor, attached to the moving body, operable to capturing video frames; and camera motion understanding device, in communication with the camera sensor, operable to determine the motion status of the moving body by determining local motion vectors from pairs of successive video frames; a classification unit operable to classify the motion status of the body as being one of a discrete number of motion status options by using mathematical autoregressive modeling of global motion parameters from the pairs of successive video frames, wherein the mathematical autoregressive models utilize a least a mean square estimator (LMS) to recursively compute the mathematical autoregressive model parameters, wherein the mathematical autoregressive modeling utilizes inputs, outputs and a white noise source, and wherein the mathematical autoregressive model is updated for every sample of the model parameter under consideration and produces a motion state classification for each model parameter by minimizing an error function.
 8. The system of claim 7, wherein the model parameters can be continuously updated over time to reflect the current status of motion of the camera.
 9. The system of claim 8, wherein the model parameters are each of: two translational components, one rotational component, two scale components, two shearing and one non-linearity component of shearing.
 10. The system of claim 8, wherein the model parameters are each of: two translational components, one rotational component, two scale components, and one shearing component.
 11. The system of claim 10, wherein the error function is time-dependent.
 12. The system of claim 11, wherein the error function associated with the LMS estimator is further adjusted for being responsive to various time-dependent conditions.
 13. A system for classifying a motion status of a moving body comprising: means for capturing video frames; means for determining the motion status of the moving body by determining local motion vectors from pairs of successive video frames; means for classifying the motion status of the moving body as being one of a discrete number of motion status options by using mathematical autoregressive modeling of global motion parameters from the pairs of successive video frames, wherein the mathematical autoregressive models utilize a least a mean square estimator (LMS) to recursively determine the mathematical autoregressive model parameters, and wherein the mathematical autoregressive modeling utilizes inputs, outputs and a white noise source, and wherein the mathematical autoregressive model is updated for every sample of the model parameter under consideration and produces a motion state classification for each model parameter by minimizing an error function.
 14. The system of claim 13, wherein the model parameters are each of: two translational components, one rotational component, two scale components, two shearing and one non-linearity component of shearing.
 15. The system of claim 14, wherein the model parameters are each of: two translational components, one rotational component, two scale components, and one shearing component.
 16. The system of claim 15, wherein the error function is time-dependent.
 17. The system of claim 16, wherein the error function associated with the LMS estimator is further adjusted for being responsive to various time-dependent conditions. 